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Deuterated molecular hydrogen in the Galactic ISM 
New observations along seven translucent sightlines 
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Abstract. We present column density measurements of the HD molecule in the interstellar gas toward 17 Galactic stars. The 
values for the seven most heavily reddened sightlines, with E(B - V) = 0.38-0.72, are derived from observations with the Far 
Ultraviolet Spectroscopic Explorer (FUSE). The other ten values are from a reanalysis of spectra obtained with Copernicus. 
In all cases, high-resolution ground-based observations of Ki and/or the CH molecule were used to constrain the gas velocity 
structure and to correct for saturation effects. Comparisons of the column densities HD, CH, CN, and Ki in these 17 sightlines 
indicate that HD is most tightly correlated with CH. Stringent lower limits to the interstellar D/H ratio, derived from the HD/2H2 
ratio, range from 3.7 x 10~ 7 to 4.3 x 10~ 6 . Our results also suggest that the HD/H2 ratio increases with the molecular fraction 
/(H 2 ) and that the interstellar D/H ratio might be obtained from HD by probing clouds with /(H 2 ) ~ 1. Finally, we note an 
apparent relationship between the molecular fractions of hydrogen and deuterium. 
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1. Introduction surements have b een performed with in the Local Interstellar 

Medium (LISM). iMoos et alJ (2002) reviewed those and re- 
It is believed that deuterium was produced in significant ported (D/H) LB = (1.52 + 0.08) xlO" 5 within the Local Bubble 
amounts only during the primordial Big Bang Nucleosynthesis. (LB) However, thi s value may not be repre sentative of the true 
Since then, deuterium has been steadily destroyed in stel- Galactic (D/HW jHebrard & ^ M^ l Eopj used a D/O survey 
lar interiors. Thus, its abundance relative to H, noted to constrain the D/H variations and, assuming a constant O/H 
D/H, is a key measurement for studies of both cosmology rati0j obtained a D/H ratio below 1.0 x 10" 5 outside the LB. 
and Galactic chemical evolution JVangioni-Flam et alj | 2000t T ndee d, recent measurements suggest a " canonical" (D/ HW 
| Coc et alj20 04). D/H in the Interstellar Medium (ISM) is char- jf it exists, likely 2 to 3 times lower dJenkins et alJ \l99%. 



acteristic of the present-day Galactic deuterium abundance. Isonneborn et allEoOol iHoones et al]l2003h . Unfortunately. ex- 
Prior to the FUSE mission, (D/H) ISM showed some dispersion tended direct invest ig at i ns over long Galactic disk sight lines 
(e.g. i Lemoine et alJ | l99Sj for a review), likely resulting from ^ difficult , due t0 saturation and blends with neighboring H 1 
poorly understood physical processes like astration, inefficient Lyman lines 

mixing, depletion onto grains, and perhaps some unidentified xhe HD/H2 ratio is an interesting alternative for D/H inves- 

systematic errors. The final resolution of these issues will have t igations along long sightlines. In the diffuse ISM, the forma- 

strong implications for our understanding of the physics of the tion of HD occurs via the i on . ne utral reaction: 
ISM, the chemical evolution of the Galaxy, and the baryonic 

density of the Universe inferred from primordial D/H. H2 + D > HD + H 

An accurate determination of (D/H) ISM is one of the main while its destruction is due to photodissociation. Because 

goals of the FUSE mission. So far, most of these mea- G f the lower abundance of deuterium compared to hydro- 

gen, self-shielding of HD becomes significant only at higher 

Send offprint requests to: S. Lacour extinction than for H2. Therefore, the transition between 
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Table 1. Target List. 



Star 


a (J 2000) 


6 (J 2000) 


V (mag) 




Distance (pc) 


Ref 


A v (mag) a 


Type 


HD 27778 


04 23 59.79 


+24 18 03.6 


6.33 


0.38 


220 


1 


1.01 


B3V 


HD 73882 


08 39 09.53 


-40 25 09.3 


7.27 


0.72 


925 


2 


2.28 


09III 


HD 110432 


12 42 50.27 


-63 03 31.0 


5.32 


0.40 


430 


3 


1.32 


B2IVpe 


HD 185418 


19 38 27.48 


+ 17 15 26.1 


7.52 


0.51 


790 


4 


2.03 


B0.5V 


HD 192639 


20 14 30.43 


+37 21 13.8 


7.11 


0.66 


1800 


5 


1.87 


08e 


HD 206267 


21 38 57.62 


+57 29 20.5 


5.62 


0.52 


615 


6 


1.37 


06e 


HD 207538 


21 47 39.79 


+59 42 01.3 


7.30 


0.64 


615 


6 


1.43 


BOV 



Ref erence s. — 1 = Simbad; 2 =lFerlet et alJ<200ct) : 3 = lvan Dishoeck & Black|jl989h : 4 = ISonnentrucker et aljiiool) : 5 jSt.-Louis & Smithl 
1199 lh : 6 = ldeZeeuwetai]<1999h . 
a Reddening and visible extinction parameters from lRachford et ail {2002 ). 



atomic deuterium and HD takes place deeper in a cloud 
than the transition between atomic and molecular hydrogen. 
Whenever all the deuterium is in molecular form, the ratio 
Af(HD)/2Af(H2) should be equal t o the elemental abundance ra- 
tio D/H (see also Wats onlll973l) . Less than a dozen observa- 
tions of deuterated molecular hydrogen have been carried out 
to date, fr om the first FUV detections with Copernicus thirty 
years ago dSpitzer et alJl973HBlack & Dalparnoll973llW atson 
19731) to infrared observations with ISO (Wright et al. 1999; 
Bert oldi et alJll999l) . Recently, the higher sensitivity of the 
FUSE satellite has allowed detection of HD in extr a-galactic 
sources dBluhm & de BoeJl200ll lAndre et alJ E()04). In addi- 
tion, more heavily reddened clouds which previously could not 
be investigated in the FUV are now within reach, as shown by 
iFerlet et all l l2000l) toward HD 73882 (A v « 2.3 mag). 

In this work, we make use of FUV data obtained with FUSE 
towards seven early-type Galactic stars (see Table Each 
sightline is "translucent", i.e., showing significant extinction 
(1 mag < Av < 5 mag; see iRachford et al1l2002t) . Because 
of the high sensiti vity of FUSE ( 10000 times more sensitive 
than Copernicus'. iMoos e t al. 2000), the data have good S /N ra- 
tios, which allow accurate measurements of equivalent widths. 
However, at the large column densities needed to have most of 
the deuterium in molecular form, most of the HD lines avail- 
able will be strongly saturated. Prior knowledge of the sight- 
line velocity structure is therefore crucial for the analysis. We 
made use of available very high resolution optic al data for 
CH - known to be a good tracer of H2 (D anks et al.l ll984; 
Mag nani et al.ll 19981) - in order to constrain the gas velocity 
structure toward our target stars. 

In the next section we describe the criteria used to select the 
targets and note some interesting aspects of each of the seven 
lines of sight. In sections 3 and 4, we present the FUV ob- 
servational data and our methods of reduction and analysis. In 
section 5, we discuss the inferred deuterium abundances. 

2. The sample 

The target selection was based on three criteria: a significant 
extinction (Ay > 1 mag), the availability of high resolution 
data for the CH molecule, and a good S /N ratio (> 10 per 
pixel) in the FUSE data. The targets of the FUSE PI 16 pro- 
gram ("Survey of H2 in Translucent Clouds"; Snow et al.) ful- 



filled the first two criteria. We retained for this survey the seven 
sightlines with the best S /N (listed in Table[Q. 

HD 27778 is located behind the outer portion of LI 506 
(Lvndsll9 62). T he che mistry of this region was investigated by 
Federma n et alJ ([1994). Their model suggests a low UV flux 
(which is likely a consequence of the filamentary structure of 
the ISM in Taurus region), as well as an average density inside 
the cloud of nearly 900 cm 3 . 

HD 73882 is very interesting for translucent cloud studies 
since it is a very bright early type star that allows us to probe 
dense material with high S IN. It is believed that this sightline is 
dominated by one o r more dense clou ds consistent with translu- 
cent cloud models ( ISnow et alJ2000T) . 

HD 110432 lies beyond the C qalsack d ark nebula. Its aver- 
age reddening (Ay — 1 .32 mag, see IRachford et alfcOOlh makes 
it intermediate between diffuse and translucent lines of sight. 

HD 185418 is a B0.5V star at 790 pc from the Sun 
dFitzpatrick & Massalll988l ll990T) . While its high reddening 
suggested the potential existence of translucent clouds along 
the sightline, a recent detailed study of the gas physical condi- 
tions showed that the sightli ne is instead compr ised of multiple 
diffuse components dSonnentrucker et al.l2003l) . 

HD 192639 is a member of the Cyg OB 1 association. As 
for HD 185418, the study of the physical conditions along the 
sightline has revealed it to be dominated by an ensemble of 
diffuse clouds (iSonnentrucker et al.l2002l) . 

HD 206267 is a quadruple system within the Cepheus 
OB association cluster Trumpler 37. Both the cluster and 
the associated H11 reg ion IC1396 have been well studied 
( Morbide lli et al.l ll997). It is believed that most of the inter- 
vening material is foreground, with a small contribution from 
IC139 6. Spectra of three of the stars in the system dPan et alJ 
2001 ) reveal substantial variations for CN on sub-parsec scales, 
but smaller variations for CH (less than 20 %). 

HD 207538 also is in the Cepheus OB2 association 
dHumphrevslll978l) . Polarization data show that this line of 
sight has a smal l Ry, but otherwise ve ry little is known about 
its ISM content dCatanzaro et all2003l) . 
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Table 2. Log of FUSE Observations. 



Star 


FUSE ID a 


Start Date 


Number of 


Exposure Time 


S/N b 








Exposures 


(ks) 


(A 1070 A) 


HD 27778 


PI 160301 


2000.10.27 


4 


9.7 


10.8 


HD 73882 


P1161301 


2000.01.24 


6 


11.9 


5.1 




P1161302 


2000.03.19 


8 


13.6 


4.6 


HD 110432 


P1161401 


2000.04.04 


5 


3.6 


28.5 


HD 185418 


PI 162301 


2000.08.10 


3 


4.4 


14.9 


HD 192639 


PI 162401 


2000.06.12 


2 


4.8 


8.1 


HD 206267 


PI 162701 


2000.07.21 


3 


4.9 


10.2 


HD 207538 


PI 162902 


1999.12.08 


4 


7.7 


6.2 




PI 162903 


2000.07.21 


10 


11.2 


7.1 



a Archival root name of targets from FUSE PI team observations. 

b Average per-pixel S/N for a 1 A region of the LIF 1A spectrum near A 1070 A. 



3. Observations & Data Analysis 

3.1. FUSE data 

The FUSE miss ion, its planning, a nd it s on-orbit performanc e 
are discussed bv lMoos et al.l J2000t) and lSahnow et alJ ( l2000t) . 

The list of the 7 targets studied in this work and the obser- 
vation information are given in Tables[2and|3 respectively. All 
data were obtained with the source centered in the 30" x 30" 
(LWRS) aperture, with total exposure times ranging from 3.6 
ks (HD 1 10432) to 25.5 ks (HD 73882). All our datasets have a 
S/N ratio per pixel between 10 and 30. All the data were pro- 
cessed with version 2.0.4 of the CalFUSE pipeline 1 . Correction 
for detector background, Doppler shift, geometrical distortion, 
astigmatism, dead pixels, and walk were applied. No correction 
was made for the fixed-pattern noise. The ID spectra w ere ex- 
tracted from the 2D sp ectra using optimal extraction JHornel 
Il986t iRobertso'r] 1 1 9861) . The ID spectra from individual ex- 
posures were cross-correlated and co-added for each detector 
segment. Equivalent width measurements were performed in- 
dependently for the spectra from each segment. Since the nom- 
inal spectral resolution is « 20 km s _I (FWHM), we binned the 
data by 4 pixels (» 7 km s -1 ) to increase the S/N ratio. All pro- 
cessed data have therefore a S/N per element greater than 20. 

3.2. Ground Based Observations 

For HD 73882 and HD 110432, high resolution CH spectra 
were obtained with the 3.6m telescope at ESO, La Silla, using 
the CES spectrograph during one run in 2001 February. The 
reduction of the data was done using a homemade IDL pack- 
age. First, we subtracted a mean bias value from each spectrum, 
then made adjustments to account for the background radiation. 
After flatfielding by means of spectra from a quartz flatfield 
lamp to remove the pixel-to-pixel sensitivity variations inherent 
to the detector, we used spectra from a Th-Ar hollow cathode 
lamp to determine the wavelength calibration. The resolution 
was estimated at 3 km s _I from the widths of the thorium lines 
in the Th-Ar exposures. 



1 http://fuse.pha.jhu.edu/analysis/analysis.html 



For the other stars, high-resolution (FWHM ~ 1.2-2.0 
kms' 1 ) spectra of K i, Na i, Ca i, Ca n, CN, CH, and CH + were 
obtained with the Kitt Peak coude feed telescope in various 
runs from 1995 to 2000. A detailed discussion of the reduction 
and analyses of these spectra - and of similar s pectra for other 
stars in the FUSE translucent cloud survey jRachford et alJ 
2002) - will be given by Welty, Snow, & Morton (in prepa- 
ration). 



4. Column Density Measurements 

Due to its dipole moment, HD is primarily detected through 
transitions from the ground rotational state 7=0. There are 20 
Lyman and 6 Wer ner HD/ = n transitions in the FUSE wave- 
length range (Dabrows ki & Herzber Jl9 76l. In Fig. [2 we plot 
the spectrum of HD 1 10432 from 1000 A to 1 1 10 A. Many ab- 
sorption lines are detected, most of them due to molecular hy- 
drogen. The upper tick marks indicate the H2 7= 1 to 4 lines and 
the thick tick marks indicate the HD 7=0 Lyman lines. Most of 
the HD transitions cannot be detected because either they are 
too close to the saturated H2 (7=0 and 7=1) absorption lines or 
they are blended with atomic lines. Therefore, only 7 HD tran- 
sitions between 959 A and 1 106 A could be used (see Table[3}- 
The wavenumbers have been determined through absorptio n 
and emission spectroscopy by Dabrowski & Herzberg lll976l) . 
Abgrall and RouefF (2004, to be submitted) have calculated the 
transition energies, the oscillator strengths, and total radiative 
lifetimes in the Lyman and Werner band systems by fo llowing a 
similar approach to the one used for H2 lAb grail et a l. 1993a b, 
2000). The rest wavelengths and oscillator strengths for the HD 
lines seen toward the seven sightlines studied here are listed 
in Table [3] Experimental wavelengths have been used when 
available; the typical accuracy on the transition wavenumbers 
is within one wavenumber. 

Two methods were used to obtain HD column densities: the 
curve of growth method (hereafter COG) and the profile fitting 
method (hereafter PF). 
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Fig. 1. FUSE spectrum of HD 1 10432 from 1000 to 1110 A. Many molecular absorption lines are detected, most of them belong- 
ing to molecular hydrogen. The upper tick marks indicate the position of H2 lines from J- 1 to 4. The thick tick marks indicate 
the position of the HD 7=0 Lyman lines. H2 lines from the ground and the first rotational states are strongly saturated, consistent 
with the high reddening in this line of sight. The gap around 1085 A is due to the gap between the two microchannel plates 
mounted in the LiF 1 A and LiF 2 A detectors. 

Table 3. Equivalent Widths of HD Lines used in this work. 



Star 


W(959) 


W(975) 


W(1011) 


W(1021) 


W(1054) 


W(1066) 


W(1105) b 




(mA) 


(mA) 


(mA) 


(mA) 


(mA) 


(mA) 


(mA) 


HD 27778 


27+11 




34±6 


34±8 


38±4 


42±4 


15±6 


HD 73882 








43±9 


45±4 


44±5 


16±4 


HD 110432 


29±3 




34±3 


32±2 


34±3 


34±3 


12±3 


HD 185418 


31±5 




58±6 


62±16 


48±6 


53±6 


21±5 


HD 192639 


48±7 


42+5 


49±4 


52±3 


53±3 




15±7 


HD 206267 


47±4 


43±4 


49±4 


53+3 


55+3 


47±3 


20±8 


HD 207538 


38±6 






45±4 


49±3 


55±4 


27±9 


Line 


L14-0R(0) 


L12-0R(0) 


L8-0R(0) 


L7-0R(0) 


L4-0R(0) 


L3-0R(0) 


L0-0R(0) 


^0 


959.82 


975.58 


1011.46 


1021.46 


1054.29 


1066.27 


1105.86 


/ 


0.0147 


0.0196 


0.0262 


0.0254 


0.0164 


0.0115 


0.000744 



a Slight blend with C 1* has been accounted for. 



4.1. Curve of Growth Method 

The equivalent widths (EqW) of the HD lines studied in this 
work are summarized in Table [5] The stellar continuum in the 
vicinity of each line was estimated using a low-order Legendre 
polynomial to fit the data. The error bars given here are l<x (see 
Sect. B"3. 



Most of the lines have similar oscillator strengths and show 
a large degree of saturation. It is therefore crucial to have ac- 
curate information on the velocity structure of the gas for these 
lines of sight. Because the integrated sightline colu mn den sities 
of CH and H? are generally very well correlate d (Danks et al. 
Il984t Ivan Dishoeck & Blacklll989t iMagnani et al.lll998l) . we 
used the CH structure to generate multi-component curves of 
growth for HD. Because Ki is generally well correlated with 
both CH and H2 and also (being heavier) exhibits smaller ther- 
mal broadeni ng, it can be used to discern even more complex 
structure JWeltv & HobbsfeOOll) . Toward HD 73882, for exam- 
ple, we used the five-component structure found for K 1 to fit 
both CH and HD. 



The multi-component curves of growth are obtained by in- 
tegrating multiple Voigt profiles, each of which corresponds to 
one gas component. Each component is defined by a fo-value, 
a radial velocity, and a column density. For each component, 
the radial velocity is taken equal to that of the CH, and the col- 
umn density is proportional to that of CH. Because HD and CH 
differ in mass (3 vs. 13 amu), we modified b as follows: 



b 2 cu + 2kT(- 



1 



1 



-) 



mno men 



(1) 



The temperature of the gas is calculated from the two lower 
rotational states of H 2 ( Rachford et alJl2002l) . Table [H lists the 
adopted component structures and the H2 gas temperature Tqi . 

Fig- 13 shows the best COG fit to the measured equivalent 
widths. The weak line at 1105.86 A is important. Its very low 
oscillator strength (almost one hundred times smaller than the 
others) makes it nearly optically thin. However, blending with 
a line due to C 1* hindered accurate measurement of its equiv- 
alent width, as reflected by the larger error bars. Nevertheless, 
because this is the only line not significantly saturated, we dou- 
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Table 4. Assumed Component Structure. 



5 



Star 




Model 


Component 


Relative 










(K) 






Strength 


km s -1 


km s~' 


kms -1 


HD 27778 


55 


CH 


1 


0.47 


14.6 


1.1 


1.3 








2 


0.53 


17.1 


2.0 


2.1 


HD 73882 


51 


Ki/CH 


1 


0.25 


20.5 


0.9 


1.1 








2 


0.45 


22.1 


0.9 


1.1 








3 


0.24 


23.7 


0.9 


1.1 








4 


0.04 


25.5 


1.1 


1.3 








5 


0.02 


28.7 


1.0 


1.2 


HD 110432 


68 


CH C 


1 


0.32 


2.9 


1.7 


1.8 








2 


0.67 


6.9 


1.3 


1.5 


HD 185418 


101 


CH 


1 


0.59 


-10.40 


1.3 


1.5 








2 


0.41 


-6.20 


2.5 


2.7 


HD 192639 


98 


CH 


1 


0.56 


-16.0 


3.7 


3.8 








2 


0.44 


-10.0 


2.1 


2.3 


HD 206267 


65 


CH 


1 


0.50 


-17.5 


2.6 


2.7 








2 


0.30 


-13.2 


1.0 


1.2 








3 


0.20 


-10.1 


2.1 


2.2 


HD 207538 


73 


CH 


1 


0.19 


-19.3 


0.8 


1.1 








2 


0.35 


-16.6 


0.8 


1.1 








3 


0.38 


-13.5 


1.2 


1.4 








4 


0.08 


-10.2 


1.2 


1.4 



11 T 0l are from the FUSE H 2 Survey iRachford et all2002h . — b See Sect.|4T] — c ICrawfordl ll995t) . 



bled its weight in the COG fits to increase its impact compared 
to that of the saturated lines. 



4.2. Profile Fitting Method 

We also used a profile fitting progr am called Owens 
jLemoineetalJl200a iHebrard et afll2002l) . developed by M. 
Lemoine at the Institut d'Astrophysique de Paris, to estimate 
the HD column densities. This program allows us both to fit 
all the HD lines simultaneously (performing a global x 2 min- 
imization) and also to add other species (e.g., Ci* and Cn** 
to deblend the 1105.86 A line). Because of the significant sat- 
uration of the other HD lines, we used the velocity distribution 
of the CH components to constrain our fit, as for the COG. 
But unlike the COG technique, profile fitting (PF) allows us 
to leave both the fo-value and relative strength as free param- 
eters for each component. This approach is more realistic, be- 
cause even though HD appears to be correlated with CH (see 
below), the correlation is not perfect. The best fits for several 
absorption lines toward HD1 10432 are shown in Fig. [3] These 
fits indicate clearly the importance of adding the Ci* 1105.73 
A line to deblend the HD 1105.86 A line; the profile fitting 
technique allows us to better disentangle the absorption from 
the two species. The derived lcr error bars are discussed in the 
following section. 

4.3. Error estimation 

We considered four basic types of errors, which are calculated 
differently depending on the method used: 



The statistical errors are assumed to be a Poissonian 
noise. Those errors (roughly the square root of the count rate) 
are computed by the CalFUSE pipeline for each pixel. For the 
COG technique, we obtained the total statistical error by sum- 
ming in quadrature the error for each integrated pixel (more in- 
forma tion can be found in Appendix A of lSembach & Savage! 
119921) . For the PF technique, this error is used by Owens to 
calculate the^ 2 . 

The background uncertainties are proportional to expo- 
sure time. The FUSE science data processing team assumed it 
to be close to 10 % of the computed background 2 . This error is 
calculated by the pipeline and added to the statistical errors. 

When using the PF technique, Owens optimizes the con- 
tinua over all the HD lines at once, so that the continuum 
placement error is included automatically in the^ 2 . The mea- 
surement of EqWs, however, requires a prior determination of 
the continuum, and in that case, the continuum error depends 
mainly upon the S /N ratio in the vicinity of the line. To esti- 
mate it, we shifted the continuum by 1 to 3 % (depending on 
the S IN ratio; see Sem bach & Savagdl 19921) . and determined 
lower and upper values for the EqW. The maximum difference 
is considered as our lcr error. 

The systematic errors are the most difficult errors to quan- 
tify. They may come from geometrical distortion, walk, dead 
pixels, point spread function (PSF) or fixed pattern noise. 
While most of these are corrected by the pipeline, small resid- 
ual effects may influence our measurements. Moreover, there 
is no way to judge their influence over a single absorption line. 
To estimate these possible systematic errors, we have assumed 
that (1) they are proportional to the quadratic sum of the pre- 

2 http://fuse.pha.jhu.edu/analysis/calfuse_wp3.html 
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Fig. 2. Curves of growth for each of the seven lines of sight. 
The strong saturation of all but one of the HD lines indicates 
that it is necessary to use accurate information on the velocity 
structure of the gas to derive yV(HD) for our lines of sight. We 
fitted our equivalent widths (EqWs) on the curves of growth 
obtained from the CH structure (see Sect. l4.1> . The unsaturated 
HD line at 1 105.86 A has a large error bar because of blends. 



viously calculated errors, and (2) the number of measurements 
performed by each method is statistically significant (i.e., they 
have a reduced^- 2 equal to 1). Under those assumptions, we 
account for systematic errors by scaling the previous error to 
have the minimum x 2 equal to the degrees of freedom. The x 1 
used for calculating the error with the PF method is implicit. 
For the COG technique, the x 1 is obtained by fitting our points 
on the multi-component curve of growth. 

There is an additional uncertainty which has the potential 
to exceed the ones already accounted for. It has to do with the 
goodness of the tracer that we used to model the HD lines. 
Because of the high S/N of our optical data, the CH component 
structures seem fairly well determined. However, it is likely 
that the HD components do not match perfectly those of CH. 
Because the degree of such differences is unknown, the result- 
ing uncertainty could not be included in the stated error bars. 
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Fig. 3. This figure shows the 2-component profile fitting of the 
HD lines toward HD 1 10432. Plotted are the HD lines present 
in the FUSE LiF 1A and LiF IB channels as well as the fit 
of the CH absorption lines from the CES data (at 3 km s ). 
Although clearly seen in the CES data, the velocity shift of 4 
km s _1 between the two CH components would have remained 
undetected using the FUSE data alone. 

Table 5. FUSE HD Column Densities. 



Star 




A'(HD)" 






FIT 


COG 


Mean 


HD 27778 


3 +2 - 9 


3 4 +34 


3 2 +3 - 2 


HD 73882 


7 4+4.9 


4 9+2.6 


5.8+^ 


HD 110432 


1 9 +0 - 5 


2 +1 

^•"-0.7 


I O+0.7 


HD 185418 


4 +L5 

^•"-0.8 


4 6 +2 2 

^•"-1.4 


4 3+1.9 


HD 192639 


| 9+O.6 
1 - z '-0.5 


1 7 +1 ' 2 

1 ' ' -0.8 


1 5 +a9 

-0.7 


HD 206267 


1 3 +L0 


9 Q+2.0 
^• y -lA 


2 1 +L5 


HD 207538 


3 5 +3 - 9 

J - -1.8 


6.6™ 


5 +5J 

j.u _ 24 



All column densities are times 10 



We note, however, that the current determination of the HD col- 
umn density allows a confirmation a posteriori of the quality of 
our tracer. 



5. Results & Discussion 

The values of yV(HD) obtained for the 7 FUSE lines of sight 
are summarized in Table [5] There is good agreement between 
the column densities determined using both the COG and the 
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PF techniques. In these particular cases, similar values are ob- 
tained by fixing only the component velocity distribution (PF) 
and by fixing in addition the relative strengths and the broad- 
ening of the CH components (COG). This agreement gives an 
indication of the reliability of our results. For the discussion 
below, we adopt the unweighted mean of the column densities 
determined by the two techniques. 

5.1. Relation between CH and HD 

The good correlation found between the column densities of 
CH and H2 in diffuse interstellar clouds is consistent with the 
production of most of the CH via a networ k of gas phase reac - 
tions in which H2 is a major component (Dan ks et alJ ll984). 
This correlation recently was extende d to more heavily red - 
dened sightlines observed with FUSE JRachford et alJl2002tl . 
The relatively tight relationship between these two species and 
the inferred chemical link both suggest that CH can be used as 
a tracer of H2. Because HD is also produced through gas phase 
reactions with H2 and because both CH and HD are subject to 
destruction through photoprocesses, one could assume that CH 
would be a good tracer for HD as well. To test this assumption, 
we plotted the column densities of H2 versus HD, and H2 and 
HD versus CH; this for all the Galactic lines of sight having 
estimates for jV(HD). T he sample includes ten sightline s ob- 
served with Copernicus ( Spitzer et alll973USnowil975l ll976, 
ll977tlWright & Morton! 1979I) (see below) and the seven FUSE 
sightlines analyzed in the present work. 

Because the A^(HD) originally given hv lSnitzeretaTllll97/!h 
were estimated assuming the HD lines (at 1054 and/or 1066 A) 
to be optically thin, those values should be taken as lower lim- 
its to the true A^(HD). We therefore have re-examined all the 
Copernicus data and have derived new estimates for A^(HD), 
using a method similar to that used for analyzing the FUSE 
spectra. The Copernicus Ul scans encompassing the HD lines 
at 1054 and/or 1066 A were obtained from the MAST archive. 
Background levels were estimated from scans of nearby sat- 
urated H2 lines obtained in the same observing programs. 
The equivalent widths measured from the normalized spec- 
tra (Table |6jl are generally in go od agreement with those in- 
ferred from the jV(H D) listed by Spitz er et all Jl973l) and with 
the values listed by Wright & Morton! il91% . The most no- 
table exception to that agreement is 10 Lac, where we sus- 
pect that the absorption near the HD line at 1054 A is due 
mostly to stellar Fe m an d/or Cr m jRogerson & Upson1ll977t 
iRogerson & Ewelllll985i) : there is no strong absorption at the 
expected position of the HD line at 1066 A. [Spitzer et al.'s 
value for A^(HD) toward 10 Lac was surprisingly high, given 
the modest E(b-v) and molecular fraction.] Component struc- 
tures derived from high-resolution spectra of CH, Ki, and/or 
Na 1 were used to model the profiles and/or equivalent widths 
of the HD lines (/il054, /I1066, and others as available). 

Our "new" Copernicus N{HD) are listed in the last col- 
umn of Table [6] and plotted with open squares in Fig. @] In 
most cases, t he ne w values are higher than those listed by 
Spit zer et alJ Jl973l) by factors of order 2; the value for ( Oph 
is about 4 times higher. On the other hand, the new value for 
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log N(H.) 




11 12 13 14 9 10 11 12 13 

log N(CN) log N(K I) 



Fig. 4. {top) N(HD) vs iV(H2) : The correlation factor close 
to 1 comfirmes a good correlatio n between the two species, as 
predicted by chemical models (Dank set alJll984l) . The slope 
is slightly greater than 1.0. — {middle) N(H 2 ) and jV(HD) vs 
JV(CH) : Tight correlation betw een N{H?) and jV(CH ) was al- 
ready observed in wide survey JRachford et aljEoo3) . and ap- 
pears in the left panel with a slope close to 1 . Af(HD) vs A^(CH) 
also shows a good relationship, but with a slope near 2.0. — 
{bottom) Af(HD) vs JV(CN) and JV(Ki) : None of the species 
exhibit a tight relationship. — Filled diamonds have A^(HD) 
derived from FUSE spectra; open squares have A^(HD) derived 
from reanalysis of Copernicus data. All the correlation factors 
and slopes are summarized in TableQ 



10 Lac is about a factor of 8 smaller, and the v alue for o Per 
is abo ut a factor of 2 smaller than that given by ISnowl d 19751 
119761) . We note that these new values could still underestimate 
the true A^(HD) if HD is actually concentrated in fewer compo- 
nents than the tracer used to model the HD equivalent widths. 

Table summarizes the correlations between several dif- 
ferent species. As shown on Fig.0] it suggests a good correla- 
tion between HD and H2, and between H2 and CH. Thus it is 
no surprise that it indicates a correlation between HD and CH 
column densities, but with a slope somewhat steeper than that 
for N{U 2 ) vs JV(CH). There are few points with JV(HD) < 10 15 
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Table 6. Copernicus HD Equivalent Widths and Column Densities. 



Star 


N(HD) a 


Inferred EqWs a 


Literature EqWs 


N(HD) 


Ref 


New EqWs 


New JV(HD) 


AA'(HD) 






1054 


1066 


1054 1066 






1054 


1066 




(dex) 


HD 21278 


14.06 


18.5 














14.2(0.1) 


+0.15 


o Per 








31.0(2.0) 


15.77 


1 






15.5(0.3) 


-0.25 


<TPer 


14.30 


32.2 




41.4(4.9) 


15.60 


2 


27.7(2.6) 




15.6(0.2) 


+ 1.30 


6 Per 


13.57 


6.0 










6.6(0.4) 




13.8(0.1) 


+0.25 


£Per 


14.15 


22.8 


16.4 








22.9(0.6) 


17.5(1.1) 


14.4(0.1) 


+0.25 


a Cam 


14.49 


49.9 


35.8 








52.0(2.0) 


48.4(2.3) 


14.9(0.1) 


+0.40 


139 Tau 


13.84 


11.2 










8.5(1.0) 




13.9(0.1) 


+0.05 


fOph 


14.23 


27.4 


19.7 


20.6(1.6) 15.9(0.8) 


14.26 


3 


18.8(1.0) 


14.6(0.6) 


14.8(0.3) 


+0.55 


59 Cyg 


13.86 


11.7 










8.8(1.1) 




13.9(0.1) 


+0.05 


10 Lac 


14.41 


41.5 


29.7 








5.5(0.7) 


[3.3(1.0)] 


13.5(0.2) 


-0.90 



References. — 1 = ISnow] <1975lll976h : 2 = ISnowl <1977l) : 3 = IWright & MortorJ ll979h 

No te. — Errors on ne w equivalent widths (EqW) are 1-sigma. 

a ISpitzer et alJ ll973l) : equivalent widths inferred from listed A'(HD) and /-values of lAllison & Dalgarnol ll969l) . assuming optically thin 
lines (see Sect. 15. II . 



Table 7. Correlation Plots. 



Quantities 






r " 


N b 


Slope c 


HD vs. H 2 






0.923 


16 


1.29±0.14 


H 2 vs. CH 






0.911 


12 


1.18+0.19 


HD vs. CH 






0.807 


13 


1.94±0.44 


HD vs. CN 






0.732 


9 


0.66+0.15 


HD vs. K I 






0.745 


14 


0.86±0.20 


HD/2H 2 vs. H 2 






0.291 


16 


0.05±0.13 


/(HD) vs. /(H 2 ) 


- (D/H)i SM 


=0.74e-5 


0.725 


16 


0.40+.10 


/(HD) vs. /(H 2 ) 


- (D/H) ISM 


= 1.52e-5 


0.725 


16 


0.19±.05 



a Linear correlation coefficient. 

b Number of points (detections only). 

c Slope and uncertainty of weighted least-squares linear fit, consid- 
ering uncertainties in both quantities. 



cm~ 2 , however, so it would be very useful to obtain more mea- 
surements of HD (and CH) at those lower column densities. We 
also investigated whether CN and K i might be useful as trac- 
ers for HD, as all three species are expected to be concentrated 
in the denser parts of diffuse clouds. In addition, Ki is gener- 
ally more readily detected than CH an d typically reveals mor e 
details of the component structure (e.g., IWeltv & Hobbsl2Q0ll) . 
Bottom panels of Fig.|4]show HD vs CN and K i. While CN and 
K i show some degree of correlation with HD, the relationships 
are not as tight as that with CH. 

5.2. The D/H ratio 

Using the Me udon PDR model (Nehme et al. to be published), 
iLe Petit et alJ J2002h studied the properties of HD in a diffuse 
cloud with n H - 500 ctrT 3 embedded in the standard inter- 
stellar radiation field. Under those conditions, HD becomes the 
reservoir of deuterium at a visual extinction of 1 magnitude, 
where the molecular fraction /(H2) of hydrogen is close to 



1. Toward our sightlines, /(H2) does not reach such a value. 
However, if we assume in each case that (1) a significant part 
of the atomic hydrogen is not associated with the observed dif- 
fuse molecular material, and (2) all the molecular material is 
in one dominant cloud, then the D/H ratio in that main cloud 
should be equal to the N(HD)/2N(H 2 ) ratio. 

Table H shows the N(UD)/2N(U 2 ) ratios for all the sight- 
lines included in this paper, i.e., our FUSE measurements and 
the values we have re-derived from the Copernicus data. All 
these values are plotted in Fig. |5J where we have added two 
refer ence values: the Lo cal Bubb le jV(D)/jV(H) ratio of 1.52 X 
1(T 5 JMoos et alJl20 02). and the ljenkins etalTll 19991) value of 
0.74 x 10~ 5 toward 6 Ori A (First de tection by [Laur ent et al. 
1979|). Several recent observations (iHebrard & MoosM /OOS ; 
Hooo es et al][2003l IWood et al. 2004) have suggested that the 
LB value may not be representative of the general ISM and 
that a "canonical" (D/H)ism (if such exists) may be closer to 
the latter value. For our sample, the average N(HD)/2N(H2) 
yields a corresponding estimate for D/H of (2.0 ± 1.1) xl0~ 6 , 
a factor 7 below the LB value and a factor 3 below the 5 Ori 
A value. Such a low value is difficult to interpret in terms of 
enhanced stellar astration or favored deuterium depletion onto 
dust grains. Indeed, given the multi-component velocity struc- 
ture of the sightlines and the many special conditions that may 
change the HD/H2 ratio from component to component, the 
most likely explanation is that we have not yet identified a pre- 
dominantly molecular cloud. 

We now turn to a more realistic view of the sightlines in 
which several molecule-containing clouds are present, with 
only some of them having all the deuterium in the molecular 
form. To some extent we can expect that the molecular frac- 
tion /(H2) will be a relevant parameter. /(H2) will increase 
with the average UV shielding and decrease with the fraction 
of smaller, more diffuse molecular clouds present. The models 
suggest that the deuterium is molecu lar only when hydro gen is 
fully in its molecular form (Fig. 1 of lLe Petit et all2002l) . This 
can be summarized over a sightline by 

/(HD)</(H 2 ), (2) 
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Table 8. Column Densities. 



Star 


E(B-Y) /(H2) '* 


iV(Hi) 


N(H 2 ) 


7V(HD) 


Ref AT(HD)/2iV(H2) yV(CH) ref 






[cm" 2 ] 

r, r+9.1 //im 


[cm- 2 ] 


[cm" 2 ] 

1+3.2 ,m 


[cm- 2 ] 

1 n <:+2.6 / ^ -> m . rv o ,™ i 



HD 73882 
HD 1 10432 
HD 185418 
HD 192639 
HD 206267 
HD 207538 



0.72 
0.40 
0.51 
0.66 
0.52 
0.64 



0.67 
0.55 
0.47 
0.32 
0.42 
0.42 



-4.8 
1 3 +05 

1 "'-0.4 

1 3+ 05 

^••-0.5 

2 0+ 08 

9 9+O.7 



(21) 
(20) 
(21) 
(21) 
(21) 
(21) 



HD 21278 0.10 0.10 b 5.5 (20) b 

oPer 0.30 0.50 8.0±2.4 (20) 

fPer 0.33 0.59 6.4±0.6 (20) 

ePer 0.09 0.21 2.5±0.5 (20) 

£Per 0.33 0.35 1.3+0.3 (21) 

a Cam 0.32 0.35 8.0±1.6 (20) 

139 Tau 0.15 0.12 8.0±1.6(20) 

£Oph 0.32 0.63 5.2±0.3(20) 

59Cyg 0.18 0.19 1.8±0.4(20) 

10 Lac 0.11 0.06 5.0±1.5(20) 



-0.8 ' 
1.3+o | (21) 

4.4* 4 (20) 

5.8 + °! (20) 



4.9*| (20) 
8.1 



0.6 

■0.7 (20) 



+0. 



3.0*> (19) 
4.0111 (20) 
4.7* 4 (20) 
33^1 (W) 
3.4!" (20) 
2.2*-» (20) 
5.4* 1 (19) 
4.4*| (20) 
2.0! » (19) 
l-7ffi (19) 



-1.7 
s+3.7 

9+0.7 



-0.6 

4 3+ 1 - 9 

1 5 +M 

2 1 +L5 
5.0i-! 



(15) 
(15) 
(15) 
(15) 
(15) 
(15) 



3+.. (13) 



1.5*! (1 4 ) 4 

3.0*» (15) 4 

4+4 (15) 4 

6! 2 (13) 4 

2.5*| (14) 4 

8+J (14) 4 

8* (13) 4 

6!j (14) 4 

8+J (13) 4 

1+2 ci 4 



3.7+ 
1.5 + 



■"• -1.4 
9 9+1.6 
^ -1.3 
9 9+0.9 

■"-0.7 
fl.7 
1.0 
+ 1.0 

-0.7 

1 4+ 10 

3 r 3 - 2 

J -i.5 



2 5 +L7 

T C+4.4 

J.8_ 20 

4 3+ 5 - 3 



9.1+ 



1+8.7 
-4.9 
3.7+ 1 - 8 



1.8 



-1.2 
3+0.9 

-0.7 

7.4!=;? 
6.8+ 7 ; 2 
2 0+ L2 
8.8^ 



(-6) 
(-6) 
(-6) 
(-6) 
(-6) 
(-6) 

(-6) 
(-6) 
(-6) 
(-7) 
(-7) 
(-6) 
(-7) 
(-7) 
(-6) 
(-7) 



3.5 ±0.4 (13) 
1.5 + 0.3(13) 
1.3 ±0.3 (13) 
2.8 ±0.5 (13) 
3.0 ±0.2 (13) 
4.3 ±0.2 (13) 



2.1 ±0.2 (13) 5 

2.2 ±0.2 (13) 5 
< 2.0 (12) 6 

1.3 ±0.3 (13) 5 
6.9 ±1.6 (12) 5 

2.5 ±0.2 (13) 5 

2.5 ±0.8 (12) 7 



Re ferences. — 1 = this pap er (FUSE data); 2 _ W elty et al., in prep.; 3 = Crawford] Il99 
iThorburn eta?] 120031) . 6 = lFederman et alJll994t) ; 7 = Thorburn, priv. comm. 

Note. — (nn) stands for 10"". 

a /(H 2 ) = 2iV(H 2 )/ (2/V(H 2 ) + /V(Hl)). — b Inferred from E(B - V). 



4 = this paper (Copernicus data); 5 



while the molecular fractions are linked together by 

/D\ = /(H 2 ) N(HD) 
\nj lSM /(HD) 2N(H 2 ) ■ 

The average value of N(HD)/2N(H2) noted above should there- 
fore only be taken as a stringent lower limit to the (D/H)ism- 
Since the two molecular fractions reflect the average effective- 
ness of the self-shielding of each species, they should be highly 
correlated, even when summing over multiple clouds. Unlike 
the top panel, the bottom panel of Fig. [5] shows an increasing 
trend for A ? (HD)/2A ? (H 2 ) as a function of /(H 2 ) - and there- 
fore a correlation between /(HD) and /(H2) - c onsistent with 
predic tions from models of the chemistry of HD jLe Petit et alJ 

Toward our highly reddened targets, atomic deuterium is 
blended by saturated hydrogen lines, and almost impossible to 
measure; direct comparison between the two molecular frac- 
tions is not feasibl e. We therefore used the two ( D/H) ratios 
noted above Jjenkins et all 19991 iMoos et all2 002 ) to examine 
/(HD) vs. /(H2) (Fig. [6ji. The increasing trend is obvious. If 
the relation between /(HD) and /(H2) can be understood, it 
could enable a wider survey of the D/H ratio via observations 
of HD. The apparent correlation between the molecular frac- 
tions should first be confirmed by observing targets in which 
both D and HD are measurable. That, however, requires tar- 
gets with low overall column density but high /(H2). Another 
possibility would be to examine individual clouds, using very 
high resolution data in the FUV (unfortunately not available at 
present). 

Further chemical modeling of diffuse clouds, including 
sums over a statistical distribution of diverse sizes and den- 



sities, could serve both to confirm the correlations and also to 
reveal other important parameters linked to the abundance of 
HD. For example, the formation rate of HD (and therefore the 
molecular fraction) is dire ctly prop ortional to the ionization of 
H and H2 by cosmic rays JFederman et all 19961 iLe Petit et alJ 
2002). It has been suggested that the observed abundances of 
H^ i n some diffuse lines of sight require a high flux of cosmic 
ravs JCecchi-Pestellini & DalgarnolEoOol iMcCall et alJl2003l 
Le Petit et alJl2004l) and that the flux of cosmic rays could be 
very inhomogeneous in the diffuse ISM (for example due to 
variation in the magnetic field). The relationship between the 
molecular fractions of H2 and HD would provide information 
on the formation rates of both molecules and on the flux of cos- 
mic rays; the variations could give some indication of the level 
of inhomogeneities in the cosmic ray flux. 

5.3. Summary 

We derived column densities of HD for 7 reddened sightlines 
using far-UV spectra obtained with the FUSE telescope. Most 
of the HD lines available in the far-UV are strongly affected by 
saturation (see the COG plotted in Fig. |3. We used high res- 
olution optical CH data to determine the velocity structure in 
these sightlines, which allowed us to correct for the that satu- 
ration. The analysis was done using both curve of growth and 
profile fitting methods. 

We combined our new results for HD with re-analyzed 
Copernicus measurements and compared the column densities 
of HD, H2, and CH. A correlation between the HD and CH col- 
umn densities is likely, but further measurements at low column 
densities are needed. 
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Fig. 5. N(HD)/2N(H 2 ) vs N(U 2 ) and /(H 2 ). The upper dotted 
line gi ves the LISM D/H value (1 .5 x 1(T 5 ) found bv lMoos et"ai] 
(2002); the lower dotted line gives the D/H value found toward 
5 Ori (0.74 x 10~ 5 ) by Jenki ns et all Jl999i) . representative of 
the lower D/H values found in several other studies. Except for 
two sigthlines, we note an increasing trend with the molecular 
fraction which does not occur as a function of the molecular 
column density. It could be an indication that iV(HD)/2iV(H2) 
did not reach the D/H value in the sightlines, showing there- 
fore the importance of examining sightlines or clouds with 
higher /(Ha). Filled diamonds have A^(HD) derived from FUSE 
spectra; open squares have N(HD) derived from reanalysis of 
Copernicus data. 



Simulations jLe Petit et all2002l) have predicted for a one- 
cloud model with Ay ~ 1 mag a iV(HD)/2iV(H2) ratio equal 
to the D/H ratio, and lower for thinner clouds. All our sight- 
lines have A v > 1, but are unfortunately composed of multiple 
clouds. We therefore obtained only stringent lower limits for 
(D/H)ism ranging from 3.7 x 10~ 7 to 4.3 x 10~ 6 . We also noted 
a correlation between the molecular fraction of H2 and the ratio 
iV(HD)/2iV(H2) which we linked to a relationship between the 
self-shielding of H2 and HD. If that relationship is confirmed, it 
would give a mean to infer (D/H)ism from observations of HD. 
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Fig. 6. /(HD) vs /(H 2 ). The left panel assumes that (D/H) ISM = 
0.74 x 10~ 5 , and the right panel assumes that (D/H)ism = 
1.52 x 10~ 5 . Both panels show a clear correlation between the 
molecular fractions. Understanding this correlation might pro- 
vide a way to deduce (D/H)ism from HD. Filled diamonds have 
Af(HD) derived from FUSE spectra; open squares have A^(HD) 
derived from reanalysis of Copernicus data (see Sect. 15. 1> . 
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